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High-statistics measurements of differential cross sections and recoil polarizations for the reaction 
7P — >■ K'^TP have been obtained using the CLAS detector at Jefferson Lab. We cover center-of- 
mass energies from 1.69 to 2.84 GeV, with an extensive coverage in the K'^ production angle. 
Independent measurements were made using the pvr" (7) and K^p{tt~ final-state topologies, 
and were found to exhibit good agreement. Our differential cross sections show good agreement 
with earlier CLAS, SAPHIR and LEPS results, while offering better statistical precision and a 300- 
MeV increase in coverage. Above « 2.5 GeV, t- and u-channel Regge scaling behavior can 
be seen at forward- and backward-angles, respectively. Our recoil polarization (Pe) measurements 
represent a substantial increase in kinematic coverage and enhanced precision over previous world 
data. At forward angles, we find that Pe is of the same order of magnitude but opposite sign as 
Pa, in agreement with the static SU{6) quark model prediction of Pe ~ —Pa- This expectation 
is violated in some mid- and backward-angle kinematic regimes, where Pe and Pa are of similar 
magnitudes but also have the same signs. In conjunction with several other meson photoproduction 
results recently published by CLAS, the present data will help constrain the partial-wave analyses 
being performed to search for missing baryon resonances. 

PACS numbers: ll.80.Cr, ll.80.Et, 13.30. Eg, 14.20.Gk, 11.55.Jy 



I. INTRODUCTION AND MOTIVATION 

Searches for missing baryon resonances in the strange 
sector have seen intense activity in recent years [IH7] , fol- 
lowing theoretical predictions that several of these miss- 
ing states have strong couplings to strange baryons [H [9] . 
For a recent review on baryons, see Ref. [10]. The 
7P — ?► K^H''^ reaction promises to play an important role 
in this regard. It is closely related to the jp — > K^A 
reaction but differs in an important aspect. Since A is an 
isoscalar, only I — ^ N* intermediate states can couple 
to JsT+A (isospin filter), while S° is an isovector, which 
allows it to couple to both I ~ ^ N* and / = | A* 
states. Coupled-channel analyses of these reactions are 
thus anticipated to be of special interest [HJ [TT] . 

The full scattering amplitude for i^+S° photoproduc- 
tion consists of eight complex production amplitudes cor- 
responding to each of the two possible spin states of the 
incoming photon, target proton and outgoing hyperon. 
Parity relations reduce the number of independent com- 
plex amplitudes to four [12] , or the number of indepen- 
dent real observables to eight, one of them being the dif- 
ferential cross section {da / d cos ) and the rest being a 
carefully chosen set of seven polarization observables [13] . 
Because of the self-analyzing nature of the hyperon (A 
and S") decays, all seven of these polarization observ- 
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ables can be measured using different target and beam 
polarization configurations. It is thus possible to com- 
pletely characterize the complex amplitude A-fp^j^+^o 
from experimental observations. With an unpolarized 
beam and an unpolarized target however, only the differ- 
ential cross sections and recoil polarizations (Ps) can be 
extracted. 

Previous high statistics measurements for 7p — >■ /C+E^ 
have been made by the CLAS [HHH], SAPHIR [17], 
LEPS [m [19], and GRAAL [20] Collaborations. The 
SAPHIR da/d cos measurements covered center-of- 
mass energies {\/s) from threshold (1.69 GeV) to 
2.4 GeV, while CLAS reached about 2.55 GeV. Agree- 
ment between the two data sets is fair, except in a few 
backward-angle bins where CLAS shows an enhancement 
around 2.2 GeV, whereas SAPHIR is more flat. The more 
recent LEPS experiment made precision measurements 
in the forward-angles for ^/s from 1.93 GeV to 2.3 GeV, 
and is in overall fair to good agreement with CLAS and 
SAPHIR. World data on polarization is consider- 

ably more sparse, however. 

In this paper, we report new measurements of 
da /d cos and the recoil polarization for — >■ X+E'^ 
using data taken at Jefferson Lab with the CLAS detec- 
tor. We have utilized the decays E° A7 and A — ptt~ , 
and have performed separate analyses for the final-state 
topologies K^pTT~ (7) and K'^p {tt~ where the final- 
state particles that are not detected are shown in paren- 
theses. Our dcr/d cos fl^jjj measurements cover center-of- 
mass (cm.) energies from near production threshold 
(1.69 GeV) to 2.84 GeV, in 10-MeV-wide bins. We also 
cover a wide angular range of —0.95 < cosO^^ < -1-0.95, 
everywhere except the extreme forward- and backward- 
angles. Our angular binning is 0.1 in cos^^J^ . For the 
recoil polarization P^, where available, we present results 
only from the K^pTT~ (7) topology, which allows us to 
preserve the spin-transfer information between the E*^ 
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and the A in the — ?► A7 decay (see Sec. 1X3). In the 



backward-angle and near-threshold bins, where statistics 
are extremely limited for the K^pn^ (7) topology, our 
Ps measurements are from the K'^p (tt^ , 7) topology in- 
stead. For the most part, these results arc in excellent 
agreement with the previous CL AS results [MJ . In a few 
cases, the polarizations reported are not consistent with 
earlier measurements. These new measurements super- 
sede the previous (lower statistics) CLAS results. 

Theoretical and phenomenological model fits to pre- 
vious world data incorporating intermediate resonances 
have typically suffered from the problem of ambiguity. 
That is, fits with different resonance contribution hy- 
potheses gave comparable values. For example, it 
has long been known that there is a "structure" at 
around 1900 MeV in both the K+A and K+E° dif- 
ferential cross sections. Mart and Bennhold |T] first 
pointed out that for K~^A, this structure could be ex- 
plained by the i?i3(1895) "missing" resonance. Janssen 
et al. 2] extended this model to X+E" and included the 
S'3i(1900) and P3i(1910) resonances as well. They found 
that Z3i3(1895) did not significantly improve the global 
fit quality for this channel and claimed that u-channel 
non-resonant processes could have a significant contri- 
bution instead. Various other groups |3l IH [6l [21] have 
incorporated several other resonances, but it is fair to 
say that none of the models are conclusive. 

Part of the problem lies in the fact that there has been 
very little data on the polarization observables, the afore- 
mentioned model fits being primarily based on cross sec- 
tion data. The problem is even more acute for the S", 
where the polarization is inherently "diluted" compared 
to the A (this point is elaborated in Sec. 1X3). High 
precision polarization measurements such as the current 
results are, therefore, much needed. With a finer binning 
in both energy and angular kinematic variables, several 
new features can be seen that were not apparent earlier. 
For instance, the SAPHIR paper [T7j noted that in ac- 
cordance with quark model predictions, their data were 
consistent with P^ « — Pa; i-e., within their measure- 
ment uncertainties, the A and TP recoil polarizations had 
comparable magnitudes and opposite signs. Our data 
show that this is not obeyed globally. The quark model 
assumes the same production mechanism for both hyper- 
ons, which can be badly broken in the resonance region if 
different resonances contribute to A and TP production. 
Given that A* states can couple only to K^YP and not 
iiT+A, this scenario cannot be ruled out. 

Secondly, the treatment of non-resonant t- and u- 
channel contributions needs to be better understood. 
The quality of previous world data in the backward- 
angles has been found to be especially poor. The present 
data demonstrate a significant rise in the cross sections 
at backward- angles, pointing strongly to w-channel con- 
tributions, and therefore lend support to the same con- 
clusion mentioned above from Ref. |2|. 

These results have the largest kinematic coverage and 
represent the most precise measurements for this reac- 



tion that are available to date. In addition to the observ- 
ables reported here, the LEPS Collaboration [T^ and the 
GRAAL Collaboration [20j have also measured photon- 
beam asymmetries for this channel, while the FROST ex- 
periment at CLAS/JLab will measure several other 
single- and double-polarization observables. The present 
work is part of a larger program within the CLAS Col- 
laboration to make precision measurements for several 
photoproduction channels '23 26J, with the goal of per- 
forming a coupled-channel partial- wave analysis (PWA). 



II. EXPERIMENTAL SETUP 

The data used in this analysis were obtained using 
real photons produced via bremsstrahlung from a 4.023- 
GeV electron beam produced by the Continuous Electron 
Beam Accelerator Facility (CEBAF) at Jefferson Lab. 
The photons were energy tagged by measuring the mo- 
menta of the recoiling electrons with a dipole magnet and 
scintillator hodoscope system [37j, resulting in a tagged 
photon energy range of 0.808 to 3.811 GeV for the cur- 
rent experiment. A separate set of scintillators was used 
to make accurate timing measurements. The photon en- 
ergy resolution was about 0.1% of the incident beam en- 
ergy and the timing resolution was 120 ps. These tagged 
photons were directed toward a 40-cm-long cylindrical 
liquid-hydrogen cryotarget inside the CEBAF Large An- 
gle Spectrometer (CLAS) detector system, which col- 
lected data events produced by scattering. Immediately 
surrounding the target cell was a "start counter" scintil- 
lator, used in the event trigger. 

Both the start counter and the CLAS detector were 
segmented into sectors with a six-fold azimuthal symme- 
try about the beam line. A non-uniform toroidal mag- 
netic field with a peak strength of 1.8 T was used to bend 
the trajectories of charged particles and a series of drift 
chambers was used for charged particle tracking. In this 
manner, CLAS could detect charged particles and recon- 
struct their momenta over a large fraction of the 47r solid 
angle. The overall momentum resolution of the detec- 
tor was ~ 0.5%. A system of 300 scintillators placed 
outside the magnetic field and drift chamber regions pro- 
vided timing information by measuring the time-of-flight 
(TOF) for each charged particle trajectory. A fast trig- 
gering and fast data-acquisition system (capable of run- 
ning at ~ 5 kHz) allowed for operating at a photon flux 
of a few times 10^ photons/s. Further details of CLAS 
can be found in Ref. [28l ■ 



III. DATA 

The specific data set that we analyze in this work was 
collected in the summer of 2004, during the CLAS "glla" 
run period. Roughly 20 billion triggers were recorded 
during this time, out of which only a small fraction cor- 
responded to K^YP events. Each event trigger required 
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a coincidence between the photon tagger Master OR 
(MOR) and the CLAS Level 1 trigger. For each charged 
particle to trigger individually, a coincidence between the 
TOF counter scintillator hit time and the start counter 
hit time for that particle was required. For the Level 
1 trigger to fire, two particles in two different sectors 
of CLAS ("two-prong" trigger) were required to trigger 
within a 150 ns coincidence time window. The final re- 
quirement was a coincidence between the tagger MOR 
and the start counter OR within a timing window of 
15 ns. Also, only the first 40 tagger counters (corre- 
sponding to the higher end of the photon energy spec- 
trum) were included in the trigger. 

During offline processing, before any physics analysis 
began, the CLAS detector sub-systems had to be cali- 
brated. This included determining the relative offsets be- 
tween the photon tagger, start counter and TOF counter 
times, as well as calibration of the drift times in the drift 
chambers and the pulse heights of the TOF scintillators. 
Energy and momentum corrections were made for indi- 
vidual particles to account for their energy and momen- 
tum losses during passage through several layers of the 
detector sub-systems. Corrections were also made to the 
incident photon energy (Ey) to account for mechanical 
sagging in the tagger hodoscope. A detailed discussion 
of the collection and calibration of this data set can be 
found in Refs. [Ml 130] • 



IV. REACTION TOPOLOGIES AND EVENT 
SELECTION 

In the reaction jp — )■ K~^YP, the E° subsequently de- 
cays as I]° — A7 100% of the time. The A further decays 
into pir^ (charged mode) with a 64% branching frac- 
tion and the rest of the time mostly into mr^ (neutral 
mode) . Since CLAS was optimized for detecting charged 
particles, we only use the charged decay mode of the A 
in this analysis. The "three-track" topology was then 
defined as jp — K^pTT~ (7) , where all three final state 
charged particles were detected and the outgoing photon 
was reconstructed from the missing momentum using a 1- 
constraint kinematic fit to zero total missing mass. This 
topology had a number of benefits. For example, the 
A decay vertex could be reconstructed from proton and 
7r~ tracking information, leading to better energy loss 
corrections for these particles. Although the detection 
of the 7r~ helped in the overall event reconstruction for 
the three-track topology, it also led to a reduction in the 
acceptance. Negatively charged particles (like n~) were 
bent inwards toward the beam line, where the CLAS de- 
tector has its lowest acceptance. 

To circumvent the above problem, we also examined 
the "two-track" topology, defined as 7p — >■ if+p(7r~,7), 
where only the if + and the proton were detected. With 
one less particle being detected, the acceptance was 
higher than in the three-track case, especially for the 
lower-energy and backward-angle regions. The two-track 



Description 


Topology 


K~''pn (7) 


K^p(-K ,7) 


Confidence level cut 


/ 




K removal cut 


/ 




Timing cuts 


/ 


/ 


Total MM cut 




/ 


Fiducial cuts 


/ 


/ 



TABLE I: List of cuts applied to the two topologies in this 
analysis. The confidence level and K'^A removal cuts used 
kinematic fitting that required the n~ to be detected and 
could be applied to the three-track topology only. The total 
MM cut came from an additional constraint on MM{K^p) 
for the two-track data set. The fiducial volume cuts were 
applied to both topologies. 

topology thus allowed for differential cross section mea- 
surements to be made from energies close to the K^TP 
production threshold and to cover almost the entire range 
in 9^^^ . Data samples for both topologies were then 
binned in 10-MeV-wide ^/s bins for further analysis. 

Each event trigger recorded by CLAS consisted of one 
or more tagged photons. To begin the event selection 
process, at least two positively charged tracks were re- 
quired to have been detected, as possible proton and K'^ 
candidates. The three-track topology required an extra 
negatively charged particle track, which was assumed to 
be a 7r~. To minimize bias, all possible photon-particle 
combinations allowed by the given topology were taken 
to be candidate signal events. Event candidates with in- 
correctly assigned photons or particle hypotheses were re- 
moved by subsequent cuts. In the following sub-sections, 
we describe each of these event selection cuts one by one, 
referring the interested reader to Ref. ^29j for further de- 
tails. Since the two topologies followed significantly dif- 
ferent analysis chains, to avoid confusion, we list the var- 
ious cuts as applicable to each of the two- and three-track 
topologies in Table |l] 

A. Confidence level cut 

Each event in the three-track data set was kinemat- 
ically fit to an overall zero missing mass hypothesis for 
the undetected outgoing photon. This was a 1-constraint 
(1-C) kinematic fit. For every event recorded by CLAS, 
both the combinations "if"*" : p : tt^" and "p : K'^ : 7r~" 
were treated as independent event hypotheses, where the 
three particle assignments corresponded to the three de- 
tected charged particles, two positively charged and one 
negatively charged. The kinematic fitter adjusted the 
momenta of each individual detected particle, while con- 
straining the total missing mass to be zero. The shifts 
in the momenta, combined with the known detector res- 
olution within the current experiment, gave a confidence 
level for the event to be the desired reaction. For a prop- 
erly tuned kinematic fitter, background events have low 
confidence levels, while real signal events populate the 
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A tof (ns) 



FIG. 1: (Color online) Confidence levels from a kinematic 
fit to the 7p — >■ K'^pTT~{'y) reaction topology. The distribu- 
tion was fairly flat above ~ 10%, as expected for real signal 
events. Background events mostly occupied the region about 
zero. These were removed by placing a loose 1% cut on the 
confidence level. 

confidence level with a flat distribution. The confidence 
level distribution is shown in Fig. [l] The peak near zero 
came from background events, which were removed by re- 
jecting any event hypothesis with confidence level < 1%. 
Above ~ 10%, the distribution was fairly flat, as expected 
for real signal events. 

B. Timing cuts 

Track reconstruction through the different CLAS de- 
tector segments yielded both the momentum p and the 
path length I from the reaction vertex to the TOF scin- 
tillator wall. The expected time-of-flight for a track hy- 
pothesized to be a particle of mass m was then given 

by 

W-^/^^- (1) 

CLAS also measured the time-of-flight, imoas, as the dif- 
ference between the tagged photon's projected arrival 
time at the reaction vertex for the given event and the 
time the given particle track hits the TOF scintillator 
wall. The difference between these two tinie-of-flight cal- 
culations gave Atof = imoas ~ ^cxp • For each track there 
was also a calculated mass rric, given by 



(a) 




calculated K* mass (GeV) 



(b) 

FIG. 2: (Color online) Timing cuts for background removal: 
(a) three-track topology, (b) two-track topology. Events ly- 
ing outside the quadrant of black curves in both figures were 
removed from further analysis. Note the logarithmic scale for 
the intensity axes. 



where /3 = l/{ct^cas)- 

Timing information in the form of Atof or was 
used to place particle identification cuts on the proton 
and tracks for both the two- and three-track topolo- 
gies. As mentioned earlier, for each pair of positively 
charged particle tracks, both "1^+ : p" and "p : K^" 
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combinations were considered and treated as indepen- 
dent event hypotheses. The cuts are shown in Fig. [2] 
where events outside the quadruplet of black curves were 
rejected. The clusters along the diagonal in Fig. [2^, were 
due to accidental coincidences with events in different 
beam bursts corresponding to the 2 ns radio-frequency 
pulses used by the CEBAF electron accelerator. In gen- 
eral, our cuts were carefully tuned to keep signal loss at 
a minimum. The only possible exception was the upper- 
left cut boundary in Fig. [2]d that was kept tighter than 
the upper-right cut-boundary to reduce the very large 
pion background. 

C. K removal cut 

Since the A and the are separated by only about 
80 MeV in mass, some iiT+A events invariably "bled-in" 
underneath the K'^YP peak. This occurred most promi- 
nently in the kinematic regions where the lab angle be- 
tween the A and TP momenta (from A'+A/if+E^ pro- 
duction, respectively) was relatively small, typically for 
high energies and forward-angle scattering. 

For the three-track topology, it was possible to effec- 
tively remove the AT"*" A contamination using kinematic 
fitting. For this, every event hypothesis was kinemati- 
cally fit to the topology 7p AT+pTr^ (nothing missing), 
which corresponded to the reaction 7p — ^ AT+A. Since 
each component of the total 4-momentum was separately 
constrained to be zero, this was a highly constrained 4-C 
fit. Events with a confidence level > 0.1% from this 4-C 
fit corresponded to AT+A background and were removed 
from further analysis. 

Fig.jsjshows the effect of this cut at high ^/s. The blue 
dotted histogram represents the unwanted A'+A events 
that were leaking in previously and has two notable fea- 
tures. The first is the long tail from the AT+A events with 
a peak at MM{K+) « 1.115 GeV. The second is that it 
shows no sign of a peak around the mass, thereby 
implying that very few good signal events were removed 
by employing this cut. The shaded histogram represents 
the AT+E" events after AT+A background removal. 

Unfortunately, this cut required the tt^ to be detected, 
and thus could not be applied to the two-track data set. 
The remnant AT+A contamination for the two-track case 
was removed during signal-background separation (see 
Sec.|V]). 

D. Total missing mass cut 

Consider the process E° A7 pn^j from the 
perspective of a 3-body reaction. If the invariant mass 
M{pTT~) is constrained to be m\, then 3-body decay kine- 
matics and the masses for the E", proton, tt^ and 7, lead 
to the bound 

0.176 GeV < M{Tr-j) < 0.251 GeV. (3) 



2.7 GeV <Vs< 2.84 GeV 




1.14 1.16 1.18 1.2 1.22 1.24 
MM(K*) (GeV) 



FIG. 3: (Color online) Effect of the A+A removal cut: for the 
three-track topology, instead of fitting to 7p — ^ K'^ptt~ (7) for 
A+E° events, we can fit to jp — >■ K^pn (nothing missing) 
for A"*" A. By rejecting events with confidence level > 0.1% 
for the latter hypothesis, we can effectively remove the A^A 
background tail (dotted histogram in blue). Only events in 
the shaded histogram were kept after this cut. 

In our reaction of interest, M(7r^7) also corresponds to 
MM{K^p). For the three-track topology, since tt~ is 
explicitly detected and the outgoing photon is recon- 
structed via kinematic fitting, the above constraint is 
satisfied nominally. 

For the two-track topology, however, non-AT+E" back- 
ground events can lie outside the bounds given by Eq. [3] 
Since MM{K'^p) also corresponds to the total missing 
mass for the two-track data set, we place the following 
additional cut for this topology 

0.16 GeV < MM{K+p) < 0.256 GeV. (4) 

These upper and lower bounds were kept slightly looser 
than the values appearing in Eq. [3] and are shown by the 
horizontal black lines in Fig. |4] 

E. Effectiveness of cuts 

The effectiveness of these cuts can be gauged by the 
percentage of "signal" events lost due to them. The 
MM( Ar+) distributions were fit to a Gaussian signal plus 
a quartic background function before and after placing 
the cuts. From this study, the loss in signal yields due 
to the cuts was estimated to be ~ 1.8% for the two- 
track and ^ 0.62% for the three-track topologies [29j . 
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We quote these as the upper hmits of the systematic un- 
certainties in our particle identification/event selection 
for this analysis. 



F. Fiducial cuts 

In addition to the above particle identification cuts, 
fiducial volume cuts were required to remove events be- 
longing to regions where our understanding of the de- 
tector performance was relatively poor. These cuts were 
motivated by differences in an empirical efficiency calcu- 
lation between the actual data and Monte Carlo which 
indicated discrepancies in the forward-angle region and 
at the boundaries of the six sectors of the CLAS detec- 
tor due to edge effects in the drift chambers. Therefore, 
events with any particle trajectory falling near the sector 
boundary regions were removed. A 0iab-dependent cut 
on cos^iab along with a hard cut at cos^iab > 0.985 re- 
moved extremely forward-going particles that coincided 
with the beam-dump direction. Localized inefficiencies 
within the CLAS detector volume such as inside the drift 
chambers were accounted for by placing trigger efficiency 
cuts on the Monte Carlo data as functions of (/)iab7 ^lab 
and \p\ for each particle track. Additional cuts were 
placed on backward-going tracks (cos^iab < —0.5). A 
minimum proton momentum cut at 375 MeV removed 
slow moving protons, whose energy losses were difficult 
to model in the detector simulation. Events with parti- 
cles corresponding to poorly-performing TOF scintillator 
counters were removed as well. 



V. SIGNAL BACKGROUND SEPARATION 

The event selection cuts were very effective in cleaning 
the data sample for both topologies. Further removal of 
background, non-iC+E" events, was affected by an event- 
based technique that sought to preserve correlations be- 
tween all independent kinematic variables [29l [3T] . The 
motivation behind this approach, as opposed to a more 
conventional sideband-subtraction method, was as fol- 
lows. 

For a reaction with multiple decays, such as in the 
present case, there are several independent kinematic 
variables (decay angles, for instance). To perform a back- 
ground subtraction, one typically bins the data in a par- 
ticular variable, such as the production angle 0^^ . This 
is because the background level can vary widely within 
the range of the kinematic variable chosen. However, 
this binning in a single variable generally does not pre- 
serve correlations present in the other independent kine- 
matic variables of interest. Therefore, one needs to bin 
the data in multiple kinematic variables, such that in 
any particular bin, the background level (both shape and 
size) remains roughly the same. Finally, the event-based 
fits using partial- wave amplitudes that we employ in the 
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FIG. 4: (Color online) In the decay E'' — >■ A7 — >■ p7r~7, the in- 
variant mass M(7r~7) is constrained to lie between 0.176 and 
0.251 GeV. For the two-track topology, M{it~^) corresponds 
to MM{K'^p), i.e., the total missing mass. Therefore, events 
lying outside the region bounded by the two horizontal black 
lines were removed from further analysis. 



Monte Carlo to calculate the acceptance for the three- 



track topology (see Sec. VI I are specifically intended to 
reproduce the correlations present in the data. There- 
fore, we have adopted a more sophisticated approach for 
background separation. 

To execute this technique for a given event, first, an 
Nc number of "closest neighbor" events were chosen in 
the phase space of all independent kinematic variables. 
Nc was typically of the order of a hundred. These Nc + 1 
events were then fitted to a Gaussian signal s(m) plus 
a background function b{m) using an event-based, un- 
binned, maximum likelihood method (the fit variable m 
being MM{K~^)). Once the functions Si{m) and bi{m) 
had been obtained from this fit for the i*'' event, the 
event was assigned a signal quality factor Qi given by: 



Q^ 



i{m^)/ {si{m.i) + bi{mi)) . 



(5) 



The Q-f actor was then used to weigh the event's con- 
tribution for all subsequent calculations. In particular, 
the signal yield in a kinematic bin with N events was 
obtained as 



N 



(6) 



Calculations were repeated using different forms of the 
background function and several different values of N^, 
without any noticeable systematic shifts in the yields. 
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The background function consisted of two parts. The 
first part was a Gaussian with a mean around the A mass 
(the width was aUowed to vary freely) - this r epre sented 
the K~^A "bleed- in" as was described in Sec. Evt. The 
second part was modeled to represent the background 
from the high-mass end. Quadratic and Gaussian vari- 
ants were tried out for this function. Trial values of Nc 
were taken as 50, 100, 150, 200 and 300. We found that 
as long as Nc was greater than w 50, the fits were stable. 
The final results presented here used Nc = 200. 

Fig. [5] shows the results from applying this method for 
the two topologies for a given ^/s bin. The background 
levels for the three-track topology varied from < 5% at 
low energies to 5-10% at higher energies, but was gener- 
ally found to be quite low. For the two-track topology, 
the background levels varied from 10-20%, the backward- 
angles having a larger percentage of background than 
the forward-angles. The total data yields after all cuts 
and background separation were ~ 4.64 million and 
~ 0.65 million for the two- and three-track topologies, 
respectively. 



(a) 



VI. DETECTOR ACCEPTANCE 
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(b) 



FIG. 5: (Color online) Signal background separation in the 
high y's regions integrated over all angles for (a) two-track 
and (b) three-track topologies. The red (shaded) histograms 
are the data weighted by Q, representing the signal, while the 
blue (dotted) histograms are the data weighted by (1 — Q), 
representing the background. Note that the three-track data 
set has a smaller background from the K'^A "tail", because 
of the additional cut described in Sec. IIVB . 



Detector efficiency was modeled using GSIM, a 
GEANT-based simulation package of the CLAS detector. 
300 million — > K~^Tt^ events were pseudo-randomly 
generated according to phase-space distributions and al- 
lowed to propagate through the simulation. An addi- 
tional momentum smearing algorithm was applied to bet- 
ter match the resolution of the Monte Carlo with the 
real data. After processing, the "raw" (i.e., original 
phase space generated) events yielded a set of "accepted" 
Monte Carlo events. The "accepted" Monte Carlo data 
then underwent the exact same series of event recon- 
struction, analysis cuts and energy-momentum correction 
steps as applied to the real data events. 

To account for the characteristics of the event trigger 
used in this experiment, two additional corrections went 
into the accepted Monte Carlo. The first of these cor- 
rections came from a trigger efficiency study using the 
7P — )■ pTT+TT" channel. This study computed the proba- 
bility that an individual particle trajectory did not fire 
the trigger, when the reaction kinematics strongly de- 
manded (via total missing mass) that the particle should 
have been present. The average effect of this correction 
was found to be 5-6%. 

The second of these corrections accounted for the 
macroscopic path length (ct ^ 7.89 cm) of the A. The 
start counter, which was included in the event trigger, 
surrounded the target cell at a distance of about 10 cm. 
In our reaction of interest, A particles would decay into 
a proton and a tt^, and these daughter particles would 
fire the trigger (only charged particles could be detected 
by the start counter). Therefore, events where the A 
decayed outside the start counter would not trigger the 
event readout. To correct for this, events in the accepted 
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Monte Carlo data set, where the A decay vertex lay out- 
side the geometrical boundary of the start counter, were 
removed from further analysis. The probability of the 
daughter proton/vr^ tracks re-entering the start counter 
region was also studied and found to negligible. On av- 
erage, the effect of this correction was about 5%. 

To form a more accurate characterization of the de- 
tector acceptance pertaining to the kinematics of the re- 
action of interest, one should use a Monte Carlo event 
generator based on a physics model, instead of a sim- 
ple phase-space generator. Typically, this is achieved in 
an iterative fashion; one starts with phase-space gener- 
ated Monte Carlo events, extracts the differential cross 
sections, fits these cross sections to a model and uses 
the model to generate new Monte Carlo events for the 
next iteration. After several such iterations, the accepted 
Monte Carlo and data distributions are expected to re- 
semble each other to a fair degree. 

However, the above procedure assumes that one has ex- 
cellent control of the signal-background separation. For a 
complicated reaction with multiple decay angles, the de- 
tector acceptance can depend on several kinematic vari- 
ables and it becomes more difficult to disentangle the ef- 
fect of the detector acceptance on signal events from that 
on the background. Our signal-background separation 
procedure, as described in the previous section, specifi- 
cally addresses this issue. By weighting every event by its 
Q- value, we are able to produce distributions of any par- 
ticular kinematic variable that include only signal events. 

For the three-track topology, we expand the scatter- 
ing amplitude A4 for the complete reaction chain 7p — 
K^T,'^ — > K^pn^j in a basis of s-channel production 
amplitudes 

n_ 

2 

M^ix,a) ^rni^'^)' (7) 

j=i P=± 

where m — {mj,mi,m-s,mf,mjf} denotes spin projec- 
tions quantized along the beam direction for the incident 
photon, target proton, intermediate E", final-state pro- 
ton and outgoing photon, respectively. The vector x rep- 
resents the set of kinematic variables that completely de- 
scribes the reaction, while the vector a denotes a set of 34 
fit parameters, estimated by a fit to the data distribution 
using the method of unbinned maximum likelihoods. The 
only assumption made here is that any distribution can 
be expanded in terms of partial waves (denoted by the 
spin-parity combination J^). Ideally, one needs to use 
a "complete" basis of such waves, but we found that 
a "large-enough" ( J-^ = i , | , . . . , ^ ) set of waves 
was sufficient to fit the data very well. The s-channel 
waves were constructed using the relativistic Rarita- 
Schwinger formalism [32^ and numerically evaluated us- 
ing the qft++ software package [33]. A full description 
of the amplitude construction and fitting procedure can 
be found in Ref. 

Based on these fit results, each accepted Monte Carlo 



event was assigned a weight /; given by, 

where we have coherently summed over the intermedi- 
ate S'^ spins. The accepted Monte Carlo weighted by 
the fits matched the data in all physically significant dis- 
tributions and correlations, as shown in Fig. [6] for the 
production angle. The detector acceptance as a function 
of the kinematic variables x was then calculated as 

/Nacc \ I A^raw \ 

^wtd(x)= (^E^J/ I E j ' (9) 

where iVj-aw and TVacc denote the number of events in the 
given kinematic bin for the raw and the accepted Monte 
Carlo data sets, respectively. 

The above procedure required knowledge of all final- 
state particle momenta. Since this was not available for 
the two-track topology, the acceptance in this case was 
calculated from the unweighted Monte Carlo intensities 
as 

^ unwtd {^X 

) = A^acc/A^raw, (10) 

where A^raw and A^acc are the same as in Eq. [9] This 
simpler expression was also used in a previous CLAS 
K^YP analysis [TS], where the effect of using a phase- 
space Monte Carlo generator, as opposed to a physics- 
model generator, was studied in detail. The conclusion 
from that study was that as long as the energy binning 
was fine enough such that the cross sections varied very 
little within each bin, the phase-space generator would 
give the correct acceptance. The previous CLAS analy- 
sis used = 25-MeV-wide bins, while our bins are even 
finer ( 10- Me V- wide in -y/s). Therefore, the above conclu- 
sion can be assumed to hold for the present case as well. 
In Sec. [Xj'V, we show that the differential cross sections 
using the two methods are in very good agreement. 

VII. NORMALIZATION 

To calculate differential cross sections, the data yields 
were normalized by the photon flux and the target factors 
as 

— ^^^^(v^,cos6lf,^.) = ( ^, r-^' . , , ) X 

y(7s,cosg5:.) 

(Acose™.)^(V^,cose^.)' ^ ' 

where At, pt, and it were the target atomic weight, den- 
sity and length, respectively, A^^i was the Avogadro con- 
stant, J-{y/s) was the photon flux incident on the target 
for the given y/s bin, A cos 0^^^ was the angular binning 
width, and y{y/s, cos 9^,^ ) and r]{y/s, cos 9^,^^ ) were the 



10 



s = 2.705 GeV 



200 



150 



CO 

I 100 



50 



■ Data 

• Acc MC (unweighted) 
A Acc MC (weighted) 



1 1 1 1 1 



I I I I I 



-0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 



cose 



FIG. 6: (Color online) Shown are the cos S^ni. distributions 
for the data, accepted Monte Carlo and accepted Monte Carlo 
weighted by the PWA fit in the ^ = 2.705 GeV bin for the 
three-track data set. Weighing by the fit results brings the 
weighted Monte Carlo distribution into excellent agreement 
with the real data. 



number of data events and the acceptance for the given 
kinematic bin, respectively. 

Photon flux normahzation for this analysis was carried 
out by measuring the rate of out-of-time electrons at the 
photon tagger, that is, hits that did not coincide with 
any event recorded by CLAS. Corrections were made to 
account for photon losses along the beam line and the 
detector dead-time. 

A separate correction to the photon flux normalization 
was required to account for the fact that only the first 
two-thirds of the photon tagger counters (1-40) went into 
the trigger. "Accidental" events corresponding to tag- 
ger counters 41-61 could trigger if a simultaneous hit oc- 
curred in the lower (1-40) counters within the same time 
window. Such "accidental" events would be triggered as 
usual and recorded by CLAS just as any other "normal" 
event. However, the photon flux calculation would not 
incorporate the associated photon corresponding to an 
invalid tagger counter. By utilizing the trigger rates in 
counters 1-40 and assuming a Poisson distribution for the 
probability of occurrence of such "accidental" events, we 
were able to correct for this feature. The boundary be- 
tween the 40*'' and 41''* counters corresponded to the en- 
ergy bin ^/s = 1.955 GeV, which had an unreliable flux 
due to this correction. As well, faulty tagger electronics 
prevented accurate electron rate measurements for pho- 
tons in the energy bins ^/s = 2.735 and 2.745 GeV [29] . 
Differential cross sections are therefore not reported at 
these three energies. However, polarization measure- 



ments do not depend on flux normalizations and are re- 
ported in these three bins. 



VIII. UNCERTAINTIES 

The statistical uncertainties for the differential cross 
sections were comprised of the uncertainty in the data 
yield and the acceptance calculation. For the «*'' event, 
the covariance matrix from the signal-background fit de- 
scribed in Sec |V| gave the uncertainty ag. in our estimate 
of the signal quality factor Qi. Summing up these uncer- 
tainties, assuming 100% correlation for events in a given 
{y^,cos9^^ ) bin, the statistical uncertainty in the data 
yield was given by 



'^data — y 




(12) 



The relative statistical uncertainty in the acceptance cal- 
culation was computed using the expression (see Sec. 5.1 
in Ref. ^) 



'1/77-1 



(13) 



Given the overall azimuthal symmetry of the detector 
about the beam direction, the data yields in each of the 
six sectors in the CLAS detector (as tagged by the sector 
in which the if + belongs) should have been statistically 
comparable after acceptance corrections. By examining 
deviations from this symmetry, we estimated the rela- 
tive systematic uncertainty in our acceptance calculation 
to be between 4 and 6%, depending on -y/s. Data col- 
lection for the present experiment occurred in bunches 
of about 10 million event triggers (called "runs"). Our 
estimated photon flux normalization uncertainty from a 
run- wise comparison of the flux- normalized K'^Tt^ yields 
was 3.2% [in]- 

In photoproduction experiments, overall normalization 
uncertainties are often estimated by comparing the total 
ttN cross sections with other world data. Since the event 
trigger for the current experiment required detection of 
at least two charged tracks, the wN channel was not avail- 
able here. However, a careful study of the cross sections 
for three different reactions {ujp, K'^A and r]p) using the 
same (present) data set in comparison with other ex- 
periments gave a flux normalization uncertainty of 7%. 
Combining this in quadrature with the uncertainty in the 
run-by-run flux-normalized yield and contributions from 
photon transmission efficiency (0.5%), live-time (3%) and 
target density and length (0.2%), we quote an overall 
normalization uncertainty of 8.3%. The last contribu- 
tion comes from the A — >■ pTr~ branching fraction (0.5%). 
A list of all the systematic uncertainties pertaining to 
da/d cos measurements for each of the two topolo- 
gies is given in Table [Hi 
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Sonrpp of TTtippTtaintv 


Topology 


K^p-K (7) 


K^p{tt ,7) 


Particle ID 


0.62% 


1.8% 


Kinematic Fitter 


3% 




Detector Acceptance 


4%-6% 


4%-6% 


Flux Normalization 


7.7% 


7.7% 


Detector Live-time 


3% 


3% 


Transmission efficiency 


0.5% 


0.5% 


Target Characteristics 


0.2% 


0.2% 


A — > p-K~ Branching Fraction 


0.5% 


0.5% 


Overall estimate 


9.7%-10.7% 


9.4%-10.4% 



TABLE II: List of systematic uncertainties for this analysis. 
The three-track topology has a lower PID uncertainty than 
the two-track topology but acquires an additional uncertainty 
from the kinematic fitting systematics. 



IX. RECOIL POLARIZATION EXTRACTION 

A. "PWA" method 

The expansion of the production amplitude using 
partial-wave analysis (PWA) techniques in Sec. VI al- 



lows for an elegant way of extracting polarization ob- 
servables for the three-track topology. We first form the 
two-component wavefunction 



I ^mim-y ) 



^-\-7ni7n-y 
A 



(14) 



out of A. 



suppressing the rest of the indices for 
the moment. The density matrix is given by 

P= ^ \i'm,7nj){i^7nimy\, (15) 

from which the expectation value of the T,^ spin in the 
direction normal to the production plane (conventionally 
denoted by y) is obtained as 



Py, = 



Tr[p] 



E 

niim. 



(16) 

The m / and m^f indices occur only in the TP decay part 
of the above amplitudes as a constant factor that cancels 
between the numerator and the denominator. The TP 
decay portion of the full amplitude in Eq.[7]can therefore 
be suppressed for extraction. 



B. "Traditional" Method 

We first define what we mean by the "helicity frame" 
of a particle. The helicity frame (HF) of any particle is 



given by an initial rotation that aligns its direction of mo- 
tion along the z-axis, followed by a subsequent boost to 
its rest frame. ^Ai/F then defined as the angle between 
the proton and the A momenta, as measured in the A he- 
licity frame, while d^^p is the angle between the normal 
to the production plane (assumed to be the y-axis) and 
the A direction, as measured in the YP helicity frame. A 
pictorial description of these two angles is given in Fig. [7] 

Fig. [7^ shows 7p — > K^YP reaction in the cm. frame, 
where the z-axis is along the beam direction and the y- 
axis is normal to the production plane. As mentioned in 
the introduction, if both the beam and target are unpo- 
larized, as in the current experiment, parity considera- 
tions imply that the induced YP polarization can only be 
along the normal to the production plane. This is shown 
by the bold arrow in red in Fig. [7^. 

To go the YP helicity frame from the cm. frame, we 
first rotate our system so that the S° momentum points 
along the z-axis and then perform a boost to the Y^ rest 
frame. Fig. jT^j shows the Y^ — >■ A7 decay in the Y^ 
helicity frame. Since the outgoing photon (shown by the 
dotted arrow) was not detected in our experiment, the 
polarization transfer from the Y^ to the A is given by 
(see Ref. [25] for a derivation) 



Pa = -Py cos 



nA 



(17) 



Note that in terms of spin structure, the — > A7 reac- 
tion is a ^ — )■ I © 1 decay, while Eq. Il7] is obtained af- 
ter averaging over the spin projections 01 the unobserved 
outgoing photon. Thus, there is a step of "dilution" in 
the accessible Y'^ spin information that occurs here. 

In the next step, we go to the A helicity frame from 
the helicity frame by making a rotation that aligns 
the z-axis with the A direction, followed by a boost to 
the A rest frame. The A — )• pTT~ decay (see Fig. [t]:) 
is a self-analyzing reaction. That is, the A polarization 
information is contained in the intensity distribution as 

I ^ {1 + aPA cos 91^ p), (18) 

where a = 0.642 ±0.013 is the A weak decay asymme- 
try |35j . Combining Eqs. 
distribution is given as 



17 and 18 the final intensity 



I ^ {1- aPy cos e^jfp cos e^^fjp). 



(19) 



Traditionally, the extraction of Py has been made using 
this intensity distribution. 

In addition to the "dilution" mentioned earlier, a fur- 
ther step of "dilution" occurs if one does not have access 
to the A momentum. It can be shown (see Appendix A 
in Ref. [TB]) that if the E^-A spin-transfer information is 
averaged over, then Eq. 19 is replaced by 



(1 - vaPycosO^^p), 



(20) 



where 



'yHF 



is the angle between the outgoing proton's 



momentum and the normal to the K^Y^ production 
plane as measured in the Y^ helicity frame and v 
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(if 7 is nof observed) 

(a) (b) (c) 



FIG. 7: (Color online) A pictorial representation of the helicity angles j? ^^'^ ^Ahf ^^'^ polarization transfer between 
the E° and the A. (a) Shows the 7p — ^ K'^Tj" reaction in the cm. frame with the y-axis as the normal to the production plane. 
The only component of the induced E" spin measurable in the current experiment is along the normal to the production plane, 
shown by the bold arrow in red. (b) Shows the — >■ A7 decay in the E" helicity frame and (c) shows the A — >■ pTv~ decay in 
the A helicity frame. See text for details. 



is a "dilution factor". Since the 7r~ from the A decay 
was not detected in the two-track topology, the A mo- 
mentum could not be reconstructed. Therefore, Eq. |20| 
applies instead of Eq. [19] for the two-track topology. 



X. RESULTS 
A. Differential cross sections 

Figs. [8] and [9] show our differential cross section results 
separately for the two topologies. The binning in ^/s is 
10 MeV and 0.1 in cos6'^j^for both cases. The three- 
track results span from 1.80 GeV to 2.84 GeV in -y/s, 
while the two-track results start closer to the threshold 
at ^ = 1.69 GeV to 2.84 GeV. The higher acceptance 
also allows greater coverage of the backward-angles for 
the two-track analysis. 

Given the widely different analysis techniques em- 
ployed for the two topologies, the agreement between the 
two results is significant. The three-track analysis made 
use of a kinematic fitter and a PWA fit-based, weighted 
acceptance calculation method, neither of which was 
available for the two-track analysis. The slight remain- 
ing differences between the two topologies are within the 
overall 10-11% systematic uncertainties. 

The consistency between results from the two topolo- 
gies allow us to quote our final differential cross sections 
as the weighted mean according to 



x(Vs, cos 9^^ ) 



cr| - 2peffO-2Cr3 -I- (t| 



where 2:2.3 are the differential cross sections and (72,3 are 
the associated statistical uncertainties for the two- and 
three-track results, respectively. Here, peff ~ 0.33 is an 
effective degree of correlation that takes into account the 
ratio between the total signal yields in the two data sets. 
The statistical uncertainty on the mean value is given by 



cr(V^, cos 6lf„j ) = 



ct| - 2pcff ^20-3 + (T| 



(22) 



The derivation of these expressions and the computation 
of pcff can be found in Ref. [25^. 

Figs. 10 and 11 show the final differential cross sec- 



(21) 



tions for the present experiment, presented at 2089 kine- 
matic points, in comparison with previously published 
high statistics measurements. The latter consists of re- 
sults from a previous CLAS experiment by Bradford et 
al. 2005 [H], a SAPHIR analysis by Glander et al. 
2004 [1^ and a set of more recent forward-angle measure- 
ments using the LEPS detector by Kohri et al. 2006 [H] . 
Overall, there is good consistency among the different 
data sets. There is a peak at ^/s ~ 1.9 GeV prominent 
over the entire angular range. In the forward-angle bins 
{cos 9 J > 0.7), there seems to be an initial dip in the 
cross section just after the 1.9 GeV peak and a subse- 
quent rise, indicative of a smaller second peak around 
-/s w 2.1 GeV. This feature was also pointed out in the 
previous CLAS analysis IS; and is present in the latest 
LEPS data 118;; as well. 

Some notable localized discrepancies also occur be- 
tween the different results. Chiefiy, this pertains to the 
"hump" at backward-angles at ^ 2.2 GeV seen in the 
previous CLAS results 15 , but not prominent in the 
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FIG. 8: (Color online) ^. (pb) vs. cos^c.m.- Differential cross section results for the two topologies in the energy range 

d cos 0^ m 

1.69 GeV < < 2.26 GeV. The centroid of each 10-MeV-wide bin is printed on the plots. Results from the two-track analysis 
are represented by red squares, and those from the three-track analysis by blue triangles. Note that we do not present results 
in the ^/s = 1.955 GeV bin (see Sec. VII I and that the y-axes are set to logarithmic scales from y^s = 2.055 GeV onwards. All 



error bars represent statistical uncertainties only. 
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FIG. 9: (Color online) ^. (pb) vs. cos^c.m.- Differential cross section results for the two topologies in the energy range 

d cos 0^ m 

2.26 GeV < < 2.84 GeV. The centroid of each 10-MeV-wide bin is printed on the plots. Results from the two-track analysis 
are represented by red squares, and those from t he th ree-track analysis by blue triangles. Note that we do not present results 
in the -^s = 2.735 and 2.745 GeV bins (see Sec. VII I and that the y-axes are set to logarithmic scales everywhere. All error 
bars represent statistical uncertainties only. 
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FIG. 10: (Color online) 



(^b) vs. ^/s in the backward-angles: final CLAS (present analysis) differential cross section 



results as the weighted average of the two topologies are in red squares. Previous CLAS results [T^ are in blue up-triangles 
while green down-triangles are results from SAPHIR [17J. All error bars represent statistical uncertainties only. 



SAPHIR data [17 . The present CLAS results, however, 
clearly confirm this structure. In fact, the SAPHIR dif- 
ferential cross sections seem to be generally lower (or 
"flatten out") towards the backward-angles, as compared 
to CLAS results, for both K+Y.^ and K+A ESI. Gener- 



ally speaking, for both hyperons, the two CLAS results 
are in very good agreement. Recent LEPS data for K'^A 
in the backward-angles [36] also shows good agreement 
with the latest CLAS K+A data mj. So it is possi- 
ble that the flattening out at backward-angles is due to 
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green down-triangles are results from SAPHIR [17] ■ The black circles represent LEPS measurements in the forward-most 
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some overall normalization issue in the SAPHIR data. 
We also note that in the two forward-most angular bins, 
the previous CLAS results |15j showed an unphysical rise 
in the differential cross sections above y^s w 2.2 GeV, 
that was attributed to systematic uncertainties in the 
acceptance (systematic uncertainties are not included in 
Figs. 10 and 11). The new results do not show this un- 



physical rise. 



B. Recoil polarizations 



As explained earlier in Sec. |IXp , going from Eq. [19| 
(three-track) to Eq. [20] (two-track) represents a second 
step of dilution in the polarization extraction procedure. 
This arises from the fact that the polarization transfer be- 
tween A and S*' in the decay E*^ — > A7 remains unknown 
for the two-track topology. The effect of this dilution 
was studied by making use of the traditional method 
of extraction for the three-track topology and comparing 
the results from using Eq.[T9]to that from Eq.|20j On the 
average, the two results agreed very well, but the scat- 
ter in the polarization data from the diluted expression 



(Eq. 20 1 was much larger. Quantitatively, by comparing 



the point-by-point ratios of the error bars, the uncertain- 
ties in the polarization extracted via Eq. [20] were found 



to be about 2.8 times larger than those via Eq. 19 



To avoid this extra step of dilution inherent in the 
two-track topology, we have chosen to present Ps results 
using the three-track data set wherever this is possible. 
This covers the energy range 1.80 GeV < < 2.84 GeV 
and the angular range -0.55 < cos9^^ < 0.95. Our 



cos flj^jjj binning is 0.1, the same as for the differential 
cross sections. However, to bolster statistics, we bin 
wider in -y/s (minimum 30-MeV-wide bins) and demand 
a minimum occupancy (Q-value weighted yields from 
Sec. |v]) of 200 events for every kinematic point at which 
a measurement is reported. 

From our discussion in Sec. |IX[ recoil polarizations for 
the three-track topology can be extracted either by the 
intensity distribution expression, or by a more sophis- 
ticated partial-wave expansion method. Fig. 12 shows 
a comparison between the two methods. The overall 
agreement is excellent, emphasizing the fact that the un- 
derlying physics extracted in both approaches was the 
same. Note that polarizations from the "PWA" method 
are constrained by Eq.[T6]to lie within the physical limits 
of ±1.0, while the "Traditional" method polarization re- 
sults are not constrained in any such fashion. We found 
that in some kinematic regions where the degree of in- 
duced S*^ polarization was sufficiently high, the "Tra- 
ditional" method sometimes gave values that were 
greater than unity. In all such cases, the "PWA" method 
gave Ps values that were close to but always smaller than 
unity, demonstrating the consistency among the two ap- 
proaches. 

The large number of parameters in the PWA fit led 
to certain difficulties in estimating the Ps uncertainties 
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FIG. 12: (Color online) Comparison of the three-track Pe 
obtained from the partial- wave expansion approach of Eq. [16] 
(red squares) and the more traditional intensity fit based ap- 
proach of Eq. 19 (blue triangles). The agreement is excellent. 



At places where Pe approaches unity, the intensity method 
polarizations tend to overshoot slightly. The uncertainties for 
the PWA method are calculated indirectly from the statistical 
spread in the data, while the traditional method uncertainties 
directly come from the covariance matrix of a fit. 



from the "PWA" method. Given the large number of 
waves employed in our PWA fits, sometimes the waves 
interfered strongly amongst each other with only a small 
contribution to the final result. However, the covariance 
matrix elements corresponding to these waves were typ- 
ically non-zero and the accumulated contribution from 
such small but non-zero covariance matrix elements gave 
rise to unphysically large estimated uncertainties. The 
uncertainties for the "PWA" method have therefore been 
obtained from the statistical spread between neighboring 
y/s measurements for a given cos bin. A similar ap- 
proach was also followed in two other previous CLAS 
analyses of the ujp [23] and if +A [35] photoproduction 
channels and found to give reasonable results. 

The estimated uncertainties from the "Traditional" 
method are simply the uncertainties obtained from un- 
binned maximum likelihood fits to the intensity distri- 
butions in Sec. |IXp . As shown in Fig. [12] the uncertain- 
ties from the two methods are comparable, demonstrat- 
ing that the likelihood fit uncertainties faithfully repre- 
sent the statistical spread in the data. In keeping with 
the internal consistency between the two methods, we 
report our final three-track recoil polarization measure- 
ments as follows. The values of the polarizations are the 
ones from the "PWA" method while the statistical un- 
certainties are the uncertainties obtained from the "Tra- 
ditional" method. We reiterate the fact here that the 
PWA expansion in Eq. [7] was specifically tuned to repre- 
sent distributions in all kinematic variables, in particular, 
the intensity distributions given in Sec. |IX|3. Therefore, 
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the two methods are completely equivalent. The advan- 
tage of the "PWA" method is that it yields values for 
Ps that are within the physical limits, while it is easier 
to estimate the statistical uncertainties using the "Tra- 
ditional" method. 

Our systematic uncertainty for Ps was estimated to be 
^ 0.03 from an examination of the systematic difference 
between the two methods. Unlike da / d cos 6^^-^ , Py, is 
bound between ±1, which is why we quote an absolute 
uncertainty here, instead of a relative percentage uncer- 
tainty. 

In the backward-angle and/or the near-threshold bins, 
where statistics are extremely limited for the three-track 
data set, the polarizations are presented from the two- 
track analysis. In general, Pj], being an asymmetry mea- 
surement, is much more sensitive to the statistics than 
the differential cross sections. This is especially true 
for the limited-statistics backward- angle/near-threshold 
bins. Therefore, we make a judicious choice when pre- 
senting the two-track polarization results, omitting mea- 
surements with unreasonably large error bars. In all, we 
present 459 individual data points. Finally, we also point 
out that while re-binning the two-track data set in y's, we 
were careful that ==1.8 GeV always lay on a bin edge. 
Since the three-track data sets extends from 1.8 GeV on- 
wards, this ensured that there was no kinematic overlap 
between results from the two topologies. 

Figs. [13] and [M] show our final measurements in com- 
parison with some earlier results from CLAS (McNabb 
et al, 2004 [IT), SAPHIR (Glander et ai, 2004 [T7]) 
and GRAAL (Lleres et ai, 2007 [20]). Previous world 
data are generally sparse and the contribution from the 
present analysis brings wide improvements in kinematic 
coverage and precision. It is noteworthy that all previ- 
ous measurements of Ps employed the diluted expres- 
sion given by Eq. |20[ Since the majority of our polar- 
ization results come from the three-track data set that 
avoids this dilution, they represent an improvement not 
just in terms of greater statistical precision, but also in 
terms of an additional physics precision. Structures that 
were hinted at by earlier measurements are now mapped 
out in much greater detail. In the three angular bins 
cos9^^ — 0.3, 0.5 and 0.7, there are a few localized dis- 
crepancies between the present results and earlier CLAS 
measurements [14 1 . We are uncertain of the precise origin 
of these discrepancies, but given that the present analysis 
incorporates a vast improvement in statistical precision, 
makes use of more sophisticated analysis tools and have 
been checked for internal consistency, the present results 
supersede the earlier measurements by CLAS. 

XI. PHYSICS DISCUSSION 

A. Differential cross sections 

The most prominent feature in our differential cross 
section results is the peak at 1900 MeV, visible over 



the entire angular range. This is not a new feature 
by itself |T5J [T7], but its interpretation in terms of s- 
channel baryon resonances has been highly controversial. 
In the introduction section, we pointed out several possi- 
ble states, e.g. Di3(1895), 53i(1900), P3i(1910), claimed 
by different groups. More recently, the Bonn-Gatchina 
group i6, 7^ has claimed that a positive-parity Pi3(1900) 
state with a two-star rating in the PDG is able to ac- 
count for features in both the differential cross section 
and polarization. There is also a possibility that this 
structure might not be a "resonance" per se, but simply 
a strangeness threshold phenomenon. The ^-A'' bound 
state mass representing the strangeness threshold is also 
around 1.9 GeV and it has been shown by Brodsky et 
al. [37H39] that at low energies, via pure gluonic ex- 
changes, QCD can give rise to an attractive Van der 
Waals force, so that cj) and N can indeed bind. It is inter- 
esting to note that recent CLAS K^Y^^ photoproduction 
data HI] also show a prominent peak at around 1.9 GeV. 
Fig. [T5| shows the differential cross sections for the three 
hyperons A, and at 9^^ = 90°. It was pointed 
out in a previous CLAS paper 41j that 7iV — > ttN scaled 
differential cross sections (s'' x da/dt) at 0^,^, — 90° also 
show a similar "rise" commencing from around 1.9 GeV. 
It is therefore possible that all these structures are con- 
nected to the same universal ss production threshold. 

The structure at y/s Ki 2.2 GeV in the backward-angles 
is also quite interesting. Recently, the LEPS Collabo- 
ration has published rjp photoproduction results at 
backward-angles that show a bump-like structure in the 
differential cross sections at around = 2.0 GeV. The 
claim has been that the bump is absent in the 77', lu and tt^ 
channels, pointing towards the conjecture that it is due 
to some resonance that couples strongly to the strange 
sector (the has a higher ss component than rf). Thus, 
it is possible that the structures seen in the 77 and S*^ 
channels are related. 

One of the long-standing issues in model calculations 
has been the contribution from background non-resonant 
processes. The additional 300-MeV coverage at the 
higher energy end provided by our results should be very 
useful in clearing up this issue. At high ^/s, where there 
are presumably few s-channel resonance contributions, 
i-channel processes are known to dominate at forward- 
angles. In the Regge description, for high s forward- 
scattering, the production amplitude should scale as 
A{s,t) ~ s~"^^\ where the Mandelstam variable t de- 
notes the exchange momentum squared and a{t) is the 
Regge trajectory of the exchanged particle (Reggeon). 
Since da/dt ~ s~^|^p, da/dt is expected to scale as 
s2(a(,ff-i)^ where Ocff is an "effective" Regge trajectory 
for multiple Reggeon exchange. Such an s~^ (signify- 
ing aoff ~ 0) scaling was already seen at lower energies 
{y/s < 2.5 GeV) in the previous CLAS analysis by Brad- 
ford et al. [15] . A very similar scaling phenomena is seen 
in the present data at the highest energies, as shown in 
Fig. [16^ . Following Guidal et al. [43 , the previous CLAS 
paper il5J also pointed out that if one assumes the t- 
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FIG. 13: (Color online) Ps vs. y/s : Recoil polarization world data in the backward- angles. CLAS (present analysis) results 
are in red squares, earlier CLAS [14) results in blue up-triangles, SAPHIR [TJ in green down-triangles, GRAAL [5^ are in 
black circles. The error bars represent the statistical uncertainties. 



channel processes to be dominated by and i\r*(892) 
exchanges, the effective Regge trajectory could be simply 
explained as aoff = ctK+ + otK* ~ for f — 0. 

The Guidal-Laget-Vanderhaeghen paper [13] also 
pointed out that a similar scaling should be observed 



at high s, u (backward-angle scattering). How- 
ever, previous world data did not have enough statistical 
precision at backward- angles to demonstrate this conclu- 
sively. Fig. [TGj j shows the u-channel scaling behavior 
in the present analysis, as demonstrated by the strongly 
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FIG. 14: (Color online) Pe vs. y^i : Recoil polarization world data in the forward-angles. CLAS (present analysis) results are 
in red squares, earlier CLAS pJJ results in blue up-triangles, SAPHIR [17] in green down-triangles, GRAAL [20^ are in black 
circles. The error bars represent the statistical uncertainties. 



collimated peak as |m| decreases. We also found that 
an scaling worked better than here. Since the 
Reggeons involved in u-channel exchanges are the A's and 
S's, aoff ~ ctA + as ~ — 1-4 at u — > 44 , which leads to 
a scaling power of ~ —4.8. Therefore, it is possible that 
the scaUng power is steeper than s~^. 



B. Recoil polarization 

The overall trend of the polarization seems to be 
that Ps is large and positive in the forward-angles and 
tends toward zero or negative values in the backward- 
angles. Many local structures are visible, especially in 



the backward-angle bins, possibly from resonance contri- 
butions, though the variations are smoother than seen 
in A'^^A |26| . In the static quark model, assuming an 
approximate SU{6) symmetry |45j . the spin-flavor con- 
figurations of the two hyperons are |A^) = \u^d^s^) and 
= \u^d^s^). Therefore, it follows that Ps « -Pa- 
A new feature that we see from the present results is 
that this prediction is explicitly broken in certain kine- 
matic regions. Fig. [T7| shows this for one region, where 
Pa [2S] and Ps are both non-zero and have the same 
sign. Similar features are visible in several other mid- 
and backward-angle bins, but not in any of the forward- 
angle bins. In other words, the SU{6) prediction is not 
observed globally. Of course, SU{6) is known to be a 
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FIG. 15: (Color online) Differential cross sections at 

Scf.m. ~ 90° for A |26) (blue squares) and (red up-triangles) 
photoproduction from proton and E~ |40| (green down- 
triangles) photoproduction from deuterium. Except for the 
extreme forward-angles, all three hyperon channels show a 
similar peak at ^/s ~ 1.9 GeV. 



broken symmetry, and it is an interesting question by it- 
self, as to why the w ~Pa prediction seems to hold 
at high -y/s and forward-angles. A possible answer may 
lie in the fact that the static quark model assumes that 
the production mechanisms for both hyperons are the 
same. This hypothesis no longer holds if A* resonances 
contribute to the if+E*' production, and SU{6) can be 
broken explicitly. 

XII. SUMMARY 

We have presented high statistics measurements of dif- 
ferential cross sections and recoil polarizations for the 
reaction jp — ^ K'^T,'^ from production threshold to 
y/s — 2.84 GeV with a wide angular coverage, using the 
CLAS detector at Jefferson Lab (electronic versions of 
the data can be found in Ref. ^46^). These new results 
significantly extend the previous ir+E° world data on 
two separate fronts. Firstly, these precision polarization 
measurements will place additional constraints on future 
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